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Spin and cyclotron energies of electrons in GaAs/Gai_ x Al x As quantum wells 
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A five-level P-p model of the band structure for GaAs-type semiconductors is used to describe the 
spin <?*-factor and the cyclotron mass m* of conduction electrons in GaAs/Gai-^ Al^As quantum 
wells in an external magnetic field parallel to the growth direction. It is demonstrated that the 
previous theory of the g*-factor in heterostructures is inadequate. Our approach is based on an 
iteration procedure of solving 14 coupled differential P-p equations. The applicability of the 
iteration procedure is verified. The final eigenenergy problem for the conduction subbands is reduced 
to two differential equations for the spin-up and spin-down states of consecutive Landau levels. It 
is shown that the bulk inversion asymmetry of III-V compounds is of importance for the spin g*- 
factor. Our theory with no adjustable parameters gives an excellent description of experimental 
data on the electron spin g*-factor in GaAs/Gao.67 AI0.33AS rectangular quantum wells for different 
well widths between 3 and 12 nm. The same theory describes very well experimental cyclotron 
masses in GaAs/Gao.74 Alo.26 As quantum wells for the well widths between 6 and 37 nm. 

PACS numbers: 73.40.Cg 73.50.Jt 73.61.Ey 
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Spin properties of electrons in semiconductor het- 
erostructures have become in recent years subject of 
intense experimental and theoretical interest because 
of their inherent scientific value as well as possible 
spintronic applications. Among numerous heterostruc- 
tures, the system GaAs/Gai-^ALjAs has won a unique 
position. First, GaAs is after silicon the most im- 
portant semiconductor material. Second, the sys- 
tem GaAs/Gai^ALrAs has well established parameters 
which allows the theorists to describe new subtle phe- 
nomena. Third, due to the advanced growth technology 
the electrons in GaAs quantum wells have very high mo- 
bilities making it possible to detect even weak effects. 
Thus it is of interest to describe precisely the spin and 
orbital electron energies in GaAs/Gai-^Al^As quantum 
wells not only for their own sake, but also as a model 
example for other heterostructures. From the theoreti- 
cal point of view, a description of magnetooptical effects 
in GaAs-type materials is not an easy task since GaAs 
is a medium gap semiconductor. As a result, its band 
structure exhibits features characteristic of narrow-gap 
materials, but it is insufficient to treat them by the mod- 
els normally used for such systems. We showed in our 
previous work on bulk GaAs and InP that the simplest 
adequate description of the band structure of these ma- 
terials is given by the so called five-level P-p model [1, 
2], Thus, the present work on the spin g*-factor of elec- 
trons in GaAs/Gai-^Al^As quantum wells was first mo- 
tivated by the fact that the theoretical description used 
in the literature for this purpose had been based on the 
three-level model [3] . Examining the problem we realized 
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that the theory [3] suffers from other deficiencies, as we 
will demonstrate below. The theory we develop for the 
spin energies is equally valid for the cyclotron (orbital) 
energies, so that we describe also the cyclotron masses, 
although here our treatment improves only slightly the 
existing approach [4]. 

An important feature of the III-V compounds is a bulk 
inversion asymmetry (BIA) of these materials. As is 
known from the fifties [5] , the BIA results in a spin split- 
ting of energies for a given direction of the wavevector 
k. It is, however, of interest to investigate how this spin 
splitting behaves in the presence of an external magnetic 
field. One can rephrase the problem by asking how the 
Zeeman spin splitting caused by the magnetic field com- 
bines with the spin splitting caused by BIA. In the bulk 
magneto-optical or magneto-transport studies one is usu- 
ally interested in the Landau levels at k z =0, since the 
density of Landau states has singularities at the vanishing 
k z . However, in a quantum well one deals with electric 
subbands for which the value of k 2 is " frozen" in the sub- 
band wavefunction along the growth direction. The BIA 
splitting is sensitive to this k 2 z value which increases with 
a decreasing well width. This in turn is reflected in the 
spin g* -factor. In the following we will be concerned with 
symmetric quantum wells, so that the Bychkov-Rashba 
spin splitting, caused by the structure inversion asymme- 
try, does not come into play [6]. 

Our calculation of the spin and cyclotron electron en- 
ergies in GaAs/Gai-^ALjAs heterostructures is based on 
the five-level P-p model (5LM). An important advantage 
of the model is that it includes the BIA mechanism of the 
spin splitting. We solve the eigenenergy problem by an 
iteration procedure which is more precise than an expan- 
sion in powers of momentum used in the literature. We 
check the precision of consecutive iteration steps. The 
results are compared with those obtained by the three- 
level P-p model (3LM) and it is shown that the latter 
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is insufficient for the GaAs-type materials. 

Our paper is organized as follows. In Section II the P-p 
theory within 5LM is formulated for the bulk semicon- 
ductor. Next the 3LM model is used to obtain results for 
the spin g*-value. It is demonstrated that the procedure 
of Ref. [3] based on the same model is inadequate. The 
iteration solution of 5LM is checked against "exact" re- 
sults for the bulk, as obtained by a numerical procedure. 
In Section III the theory for the quantum wells is worked 
out from the 5LM matrix. Results for the spin g*-factor 
and the cyclotron mass are presented, compared with ex- 
perimental data and discussed in Section IV. The paper 
is concluded by a summary. 



II. P p THEORY 

In order to discuss the previous treatments and to es- 
tablish standards of precision for various approximations 
we first consider the 3D case. The P-p theory, which is 
the kp theory generalized for the presence of an external 
magnetic field B, has the form (cf. Ref. [1]) 
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FIG. 1: Five- level model for the band structure of GaAs-type 
semiconductors. Energy gaps, spin-orbit energies, interband 
matrix elements of momentum and of the spin-orbit inter- 
action are indicated. Letters C and C mark symbolically 
far-band contributions to the effective mass and the spin g*- 
factor of conduction electrons, respectively. 
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where £ is the energy, P = p + eA is the kinetic momen- 
tum, A is the vector potential of magnetic field B, and 
o'l'l = (1/^) < u v W\ u i >■ Here er are the Pauli spin ma- 
trices, fi is the volume of the unit cell, ui are periodic am- 
plitudes of the Luttinger-Kohn functions, hb = eh/2rriQ 
is the Bohr magneton, pin are the interband matrix el- 
ements of momentum and Hf,° are those of spin-orbit 
interaction. The sum in Eq. (1) runs over all bands 
I = 1,2, ... included in the model, I' = 1,2, ... runs over 
the same bands, E^ are the band-edge energies, see be- 
low. 

In the following we shall be concerned with the five- 
level model of the band structure, as illustrated in Fig. 
1. Within this model there exist three interband matrix 
elements of momentum 
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and three matrix elements of the spin-orbit interaction 
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where the nonprimed functions are related to the Tg 
valence bands, while the primed functions are related 
to the r§ conduction bands. The band edge ener- 
gies are generally influenced by A (cf. Ref. 1) but 
in GaAs-type materials, where A is small, this feature 
may be ignored and we have: E^ = E^ = E^ = 
B< 9 ) = Gi = E 1 + Ai; E^ = E^ = E i; E^ = 

E {6) = £(12) = £(13) = Eq . E (7) = £(14) = Gf) = 

Eq + Ao- We take the zero of energy at the Tg con- 
duction band edge. Within 5LM there arc 14 basis 
functions and the P-p differential Hamiltonian for the 
envelope functions //(r) has the explicit form of Eq. 
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= P /h, P 1 = Pi/h, Q = 
.ft - 



Q/h, 



(8). 

A = £ - hw° c (a+a + ±) - p 2 z /2m . Here uj° c = eB/m 
is the cyclotron frequency with free electron mass. The 
quantities /6 ± = C[huj°'(a + a + ±) + p 2 z /2m ] ± C \i B B 
result from the far band contributions to the conduction 
band energies. The operators P± = [P x ± iP y )/\/2 are 
proportional to the raising and lowering operators for 
the harmonic oscillator functions: P + = —{h/L)a + and 
P_ = —{h/L)a, in which L = {h/eB) 1 / 2 is the magnetic 



radius. The matrices A\, A2 and D have different forms 
for various directions of magnetic field. In the follow- 
ing we shall be concerned with the case B|||[001] crystal 
direction, for which 



Ai = A 2 = 



fie p. 
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For Q=0 and k z —0 matrix (8) factorizes into two 7x7 
matrices containing spin-up and spin-down electron en- 
ergies. For Q = and k z ^ the resulting bands are 
spherical (i.e. the energies do not depend on the direction 
of B), but the spin- up and spin-down states are mixed 
by the k z terms. Still, in this case one can find solutions 
of the eigenvalue problem in form of single harmonic os- 
cillator functions. If Q 7^0, the resulting bands are non- 
spherical and the above simple solutions do not exist. 
At this point we emphasize that matrix (8) contains no 
approximations within 5LM. 

For considerations of the bulk there is no external 
potential. To describe the magnetic field we take the 
gauge A=[-By, 0, 0], corresponding to B||[001]. One 
can then look for the envelope functions in the form 
fl = exp(ik x x+ik z z)$i(y), where $/(?/) are the harmonic 
oscillator functions. 



A. Three-level model 



The three level model of rg, T v 7 , T$ levels (sec Fig. 
1) is not adequate for the description of the conduction 
band in GaAs - type materials (see [1, 2]), but it can be 
solved exactly and we use it as a starting point of our 
considerations as well as for a discussion of the existing 
work on the g* -values. 

Once the higher conduction levels are omitted, the cou- 
plings Pi, Q and A do not come into play (see Fig. 1), 
the conduction band is spherical and the resulting 8x8 
Hamiltonian may be solved in terms of eight harmonic 
oscillator functions. This was first done by Bowers and 
Yafet [7], see also [8, 9]. In our more general formulation 
(8), the 3LM corresponds to considering the columns and 
rows 4, 5, 6, 7, 11, 12, 13, 14. The resulting energy for 
the spin-up state of the n-th Landau level is given by 
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while the energy for the spin-down state of the n-th Lan- 
dau level is given by 
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spin-up and spin-down energies one can find the cyclotron 
effective mass and the spin g*-value exactly within 3LM. 
We define the cyclotron mass m* and the spin g*-value 
in the standard way 



c± _r±_ 



heB 



ml 



£- 



(14) 



(15) 



The dependence of m* and g* on the electron energy 
was established in the early days, see Refs [11, 12]. Know- 
ing the solutions of Eqs. (11), (12) and fixing magnetic 
field intensity B one can determine the dependence of m* 
and g* on k 2 exactly within 3LM. 

Since the conduction band in GaAs is only weakly non- 
parabolic, we can expand the exact energies using Eqs. 
(11) and (12) in the limit £/E « 1. The zero-order 
terms give uiq and <?q values at the conduction band edge. 
The first-order terms give the first nonparabolic approx- 
imation. One can alternatively say that the zero-order 
terms are proportional to Pq while the first-order terms 
are proportional to Pq. Using Eqs. (11), (12) and the 
definitions (14) and (15) we obtain after some manipu- 
lation (we keep the free electron terms and the far-band 
contributions only in the conduction band) 
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G\ = Go - \» B B + hj° c (n -\)-e , 

G% = Go + \^ B B + hoj° c (n +\)-E . (13) 

The first terms in Eqs. (11) and (12) represent the orbital 
parts, the second are related to the spin parts (they differ 
in sign for spin-up and spin-down states) and the third 
give k z parts for the motion along the magnetic field. 
Equations (11) and (12) represent fourth-order polyno- 
mials for the unknown energies. If we neglect the free- 
electron and far-band contributions, Eqs. (11) and (12) 
reduce to cubic equations for the energies. An influence 
of free electron terms on the description of spin g*-factors 
in selected materials was discussed by Singh et al [10]. 

In investigations of bulk semiconductors one is usually 
interested in the simplified case of k z = since the sin- 
gular density of states for the Landau levels corresponds 
to the vanishing k z . However, our final aim in this work 
is to investigate quantum wells, which corresponds to the 
situation of k z ^ 0. Solving Eqs. (11) and (12) for the 
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where D nkz = huj° c (n + 1/2) + h 2 k 2 /2mo 1 and 

are the band-edge values of the effective mass and of the 
g*-factor according to 3LM. We define Ep = 2moP 2 /h 2 . 
To calculate the numerical values we take the following 
band parameters for GaAs: Ep = 27.865cV, E = -1.519 
eV, G = -1.86 eV, C = -3.070, C = -0.102. This results 
in m* = 0.0660 m [2] and gfi = -0.44 [13]. 

We would like to plot the g* value as a function of 
energy in view of the applications to quantum wells 
(QW). To this end we calculate the theoretical en- 
ergy £ for the ground electric subband in a rectangular 
Gao.67Alo.33As/GaAs/Gao.67Alo.33As QW as a function 
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FIG. 2: Theoretical spin (j*-factor of conduction electrons in 
bulk GaAs versus the wave vector k z o parallel to the mag- 
netic field B, or the corresponding well width d (see text). 
The solid line (EX) shows exact results from the three-level 
model, the dashed line (ap) represents the first nonparabolic 
approximation to 3LM, the dashed-dotted line (IK) illustrates 
the results of Ref. [3] also based on 3LM (see text). 



of the well width d. Next we define the value of k z o by 
the equality: £{d) = H 2 k 2 /2rriQ, where = 0.0660 mo- 
This value of k z0 = {2m>l£ /h 2 ) 1 / 2 is introduced to the 
above relations for g* . When describing the final results 
for the <7*-value and the cyclotron mass in QW we will 
use the same range of well widths. 

In Fig. 2 we plot the spin g*-value of the conduction 
electrons in bulk GaAs as a function of d (or k z0 ). The 
solid line indicates the exact values of g* within 3LM, as 
obtained from Eqs. (11) and (12) for B = 2 T. The free- 
electron and distant band contributions are kept only for 
the conduction band. The dashed line indicates values 
obtained from Eq. (17), resulting from the expansion 
for weak nonparabolicity. As can be expected, for big- 
ger widths d (or smaller values of k z o) the approximated 
values approach the exact ones. Finally, the dashed line 
marked with IK indicates values obtained from Eq. (9) 
of Ref. [3] with the use of the same material parameters. 
It can be seen that Ref. [3] predicts a much stronger 
dependence of g* on k z o than both our curves and, in 
particular, our dashed curve which has been calculated 
in the same nonparabolic approximation as that used in 
Ref. [3] . The reason of this discrepancy is that the formu- 
las given in Ref. [3] correspond to our Eq. (17), with our 
2to /too replaced by — 2E P (l/G + 2/E )/3 and our g$ 
replaced by ~2Ep(l/Go — l/£o)/3. In case of the mass 
the approximation used in Ref. [3] is not too bad since it 



misses only the term of 1+C [cf. our Eq. (18)]. In case 
of the g*-value, however, it is seriously inadequate since 
it misses the additive term of +2+2C, which is essential 
for the band-edge value of g$ in GaAs, see Eq. (19). As a 
result, the second term in the square bracket of Eq. (17) 
appears not with the coefficient g$ = -0.44, as it should, 
but with the coefficient of about -2.24. This gives a much 
too strong dependence of g* on k zQ , as illustrated in Fig. 
2. 

The dependence of g* on d (or k z o), as given by the 
exact solution of 3LM (curve EX in Fig. 2), does not de- 
scribe experimental data in GaAs/Gao.67Alo.33As quan- 
tum wells. In particular, it docs not exhibit the change 
of sign to positive g*-factors for small values of d. On the 
other hand, we know that the inclusion of higher conduc- 
tion Fy, Tg bands in the P-p description (see Fig. 1) re- 
sults in stronger conduction band's nonparabolicity and, 
in particular, it gives the electron g*-values that change 
sign to positive at higher electron energies [1] . Thus it is 
necessary to include the Ty, T| levels in the refined de- 
scription of g* in GaAs. This is equivalent to using the 
full matrix (8) for the determination of electron energies. 
In view of the applications to quantum wells the k z terms 
must be retained. We pursue this program below. 



B. Five-level model 

In this section we calculate the spin g*-valuc of conduc- 
tion electrons using the five-level P-p model. This model 
is adequate for the GaAs-type materials. Once the two 
upper conduction levels are included, the matrix elements 
Pi, Q and A come into play, see Fig. 1. As a result, two 
qualitatively new features appear. First, the appearance 
of the matrix element Q makes the conduction band non- 
spherical. Second, the matrix element Q does not vanish 
if the crystal is characterized by the bulk inversion asym- 
metry (BIA) . This results in the spin splitting of conduc- 
tion energies for a given direction of the wave vector k 
(the Dresselhaus mechanism [5]). We develop a workable 
approach to the calculation of g*-factor in quantum wells 
based of 5LM using an iterative expansion in powers of 
the interband matrix elements P , Pi and Q. Clearly, 
it is necessary to check whether such an expansion is 
adequate. In other words, we should make sure that the 
expansion procedure gives the c/*-values comparable with 
those obtained from the complete matrix (8) . To this end 
we first solve the eigenenergy problem given by the P-p 
Hamiltonian (8) by an alternative method, first used by 
Evtuhov for the nonspherical valence bands of Ge [14]. 
As already mentioned, for spherical energy bands the so- 
lutions of the eigenenergy problem in the presence of a 
magnetic field are given in terms of simple harmonic os- 
cillator functions. For not too strong deviations from 
band's sphericity we look for the envelope functions in 
the form of sums of harmonic oscillator functions with 
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unknown coefficients 

f l (r)=exp(ik z z)J2c l Jm> . (20) 

m=0 

Since our Hamiltonian (8) contains raising and lowering 
operators for the harmonic oscillators, it is possible to 
perform the prescribed operations. Then we multiply 
the obtained eigencnergy equations on the left by con- 
secutive harmonic oscillator functions and integrate the 
scalar products using the ortonormality relations. This 
results in equations for the coefficients and the condition 
for nontrivial solutions gives the eigenenergies. In other 
words, the Evtuhov method is a standard way of trans- 
forming a differential eigenvalue problem into an alge- 



braic problem by taking for the complete set of functions 
the harmonic oscillator functions. 

Not going into details of the procedure we will give 
a general scheme for the resulting eigenenergy matrix. 
The elementary irreducible block is a 7x7 matrix, one 
containing the spin-up state and another the spin-down 
state. If Q and k z terms are neglected, the big matrix 
factorizes into the 7x7 matrices for the spin-up and spin- 
down energies for consecutive Landau levels. The nondi- 
agonal Q and k z terms couple the elementary 7x7 ma- 
trices. We find that in order to obtain precise energies 
for the lowest Landau levels we need to go to 70 x 70 ma- 
trices. The schematic form of such a matrix is given in 
Eq. (21). 
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We use a symbolic notation, for example N, + stands 
for the elementary 7x7 matrix containing the conduction 
Landau level n with spin up. The symbols Pk z and Qk z 
stand for 7x7 matrices containing Pok z , P\k z and Qk z 
terms, respectively, while QV~B stands for 7x7 matrices 
containing Qy/B terms. 

To calculate the g*-value we will be interested in the 
lowest conduction levels 0^ . For N=0, matrix (21) gives 
the energy of the conduction level + (second row and 
column) perturbed by interactions with other levels. For 
N=l we obtain the energy of the 0~ level (first raw and 
column). Fixing the values of B and k z we obtain from 
matrix (21) the energies which we consider to be "exact" , 
so we use them as standards in the estimation of following 
expansions. 

Next we turn back to the set of 14 equations for the 
envelope functions ft, as given by the initial differential 
matrix (8). We find the envelope functions by substi- 
tution using an iteration procedure. In the first step we 
put Q=0 and express twelve envelope functions by ft\ and 
/n, corresponding to the spin- up and spin-down Y% con- 
duction states. Next the Q terms are restored and the 
zero-order functions are put back to the complete equa- 
tions. The twelve first-order ft functions are expressed by 
fi and /n. In this approximation, linear in Q terms, the 
BIA (Dresselhaus) spin splitting is included. Finally, we 



TABLE I: Band parameters of Ga.i- X ALAs alloys for different 
chemical compositions x, as used in the calculations, see [13, 
15-18]. Energies E , Go, Ei, Gi are defined in Fig. 1, C 
and C are far-band contributions to the band-edge values 
of mo/mS and go, respectively, see Eqs. (18) and (19). The 
interband matrix elements of momentum and of the spin-orbit 
interaction are taken to be independent of x: _Ep =27.865 eV, 
£ Pl =2.361 eV, £ Q =15.563 eV, A=-0.061 eV. 





GaAs 


Gao.74Alo.26As 


Gao.67Alo.33As 


Gi(eV) 


3.140 


3.395 


3.463 


Ei(eV) 


2.969 


3.221 


3.289 


Eo(eV) 


-1.519 


-1.888 


-1.992 


Go(eV) 


-1.860 


-2.199 


-2.297 


C 


-2.297 


-1.904 


-2.248 


c 


-0.025 


-0.043 


-0.056 


Vs(eV) 




0.208 


0.264 


mS/mo 


0.066 


0.0803 


0.0875 


go 


-0.44 


4-0.40 


+0.54 



return again to the initial set and determine the twelve 
functions ft by f& and /n including the Q 2 terms. The 
resulting equations for f± and /n represent the effective 
Hamiltonian. It is obtained in the form of a 2x2 differ- 
ential matrix since the bulk inversion asymmetry mixes 
the two spin states. 
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FIG. 3: Theoretical spin g*-factor of conducting electrons in 
bulk GaAs versus the wave vector k z o, or the corresponding 
well width d (see text), as calculated from the five- level model. 
I - results for Q = (first iteration step) , I+II - bulk inversion 
asymmetry included (linear Q terms, second iteration step), 
I+II+III - quadratic Q terms included (third iteration step). 
The line marked EX shows "exact" results calculated numer- 
ically for the bulk by the Evtuhov method. 



We are now in a position to calculate electron energies 
and the resulting conduction g*-factor for given values of 
B and k z (still for the bulk). In the obtained formulas for 
energies there appear k z , k z and k z terms. We checked 
that the k z and k z terms have negligible influence on the 
final energies. 

In Fig. 3 we plot the calculated <?*-values as functions 
of d (or k z o) for the three consecutive approximations. To 
determine k z0 for a given width d we use the same proce- 
dure as described before (cf. Fig. 2). Our three approx- 
imations for the g*-value correspond to the three itera- 
tion steps. Curve I corresponds to Q=0, the correspond- 
ing formulas contain Pq and P 2 type of terms. Curve 
I+II includes the previous approximation plus the terms 
PoPiQ, APq Q, AP 2 Q, accounting for the spin splitting 
due to BIA. Curve I+II+III includes the two previous 
approximations plus the terms PqQ 2 , P 2 Q 2 , AP PiQ 2 
(the last terms are very small and may be neglected) . In 
the same Fig. 3 we plot the g* value obtained from the 
"exact" numerical procedure based on matrix (21). 

A discussion of the results shown in Fig. 3 is in or- 
der. First, curve I gives the results very similar to those 



marked by EX in Fig. 2. Although the calculation I in- 
cludes the higher conduction levels, as compared to 3LM 
used for EX, it turns out that the results of 3LM and 
5LM with Q=0 are very similar if one adjusts C and C 
constants to have the same band-edge values of and 
<7o- Second, sizable difference between the calculations I 
and I+II may appear unexpected. Compared to I, the 
calculation I+II includes the Dresselhaus spin splitting 
due to BIA, which is usually rather small. Indeed, for 
relatively large d (small k z0 ) the corrections due to BIA 
are small. It is at small widths d (large k z o) that they 
become appreciable. One should keep in mind that BIA 
works directly on the spin splitting, i.e. on the g* value. 
More generally, the BIA spin splitting has somewhat par- 
ticular properties since, while it is affected by a magnetic 
field, its primary origin is not a magnetic field. Finally, 
it follows from Fig. 3 that the curve I+II+III almost 
coincides with the "exact" result. This means that the 
third step in the iteration procedure gives enough preci- 
sion in the calculation of the electron g*-value. We use 
this information in the following treatment of the spin 
splitting in quantum wells. 



III. ^'-FACTOR AND CYCLOTRON MASS IN 
QUANTUM WELLS 

We apply the above results to calculation of the spin 
g*-factor for conduction electrons in QWs grown along 
the z direction. We use 5LM and consider both the 
growth direction and a magnetic field parallel to the [001] 
crystal axis. Matrix (8) must now be completed by the 
potential V(z) in all diagonal terms. The calculation is 
carried out by iterating solutions of Eq. (8), as described 
above for the bulk. The difference is that now we keep the 
operator p = (h/i)d/dz in the differential form since the 
envelope functions in QWs are unknown and they must 
be determined by solving the cigenencrgy problem. As a 
consequence, it is necessary in the iteration procedure to 
observe the order of ^-dependent terms. In our treatment 
we keep the frcc-clcctron and far-band contributions only 
in the conduction band. While a more general treatment 
is possible, it makes the final formulas very lengthy while 
the precision is not markedly higher. 

After a considerable manipulation the effective Hamil- 
tonian for the spin-up and spin-down electron states in 
the conduction band is obtained in the form 



H 



A+ K 



(22) 



where 
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A+ = v(z) 



h 2 d 1 d Pi + Pi n B B 
Ydzm*{£^z)dz + 2m*{£,z) + 2 



— [S Po (X 5 P_P + P_P + + ^ 6 p2p2 + Kl p2_pT ) + E Pl (K 8 P_P + P_P+ + KvP\Pl + K w P 2 _Pl)] I , (23) 



A- Viz) 1 * + P ' + ^ - ^1 Q *(£ z) 

~ [ ) 2 dzm*(£,z)dz + 2m*(£,z) 2 j 



+ 



ft 2 V En ( 



2m J m Xdz l^(K,P + P- + K 2 P.P +) + E Pl (K 3 P + P. + K i P.P +)] - + 



2 



[E Po {K 5 P+P-P + P- + K 6 P 2 _Pl + K 7 PlP 2 _) + E Pl (K 8 P+P-P + P- + K 9 P 2 Pl + K w PlP 2 )] J 



The effective mass is 
too 



to|(£,z) 1 + C 3 



-Bo Go 



and the g*-factor is 



9l (£,z) = 2 + 2C' + 2 - 



E\Gq EqG\ 



Ep , i _ ±) +E j ± iv^f^^ 



-E-o Go 



Gi i?i 



Pi Go EqGi 



(24) 



(25) 



(26) 



where ^ 1 

A 4 = — 



= Si - 5 + , (27) 



G, = G 4 -£ + V(z) . (28) 



Gj \G £o/ Pi-Eq \Gi i?i 



1/1 1 N ' 



K 5 

3Gi \Eq Go 



1 / 1 



We use the notation E P(i = 2H 2 Pq /too, = 2h 2 P 2 /m ^ 
and E Q = 2h 2 Q 2 /m . Further 6 2£g V-Bi G 1 

1/2 1 

2 IgTg§ + e~¥ 

E 2 \Ei Gj ' GiGq \E G J K s = — [jr ~ 



1 / 26 5 5 \ 
K 1 = — — + — + — , 1 

3Gi V Go-Bo G 2 E 2 ) K 7 

1/8 1 \ 1 / 2 1 \ 1/1 ix2 
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K w = 



2 \E E 2 G Q G\ 



(29) 



We emphasize that the mass m} and the g-value g*[ 
defined in Eqs. (25) and (26) do not represent the final 
cyclotron mass and g-value in a quantum well but only 
the first iterative approximations to these quantities, as 
obtained from matrix (8) by putting Q = 0. 

The nondiagonal component in Eq. (22), related to 
the bulk inversion asymmetry, is 



K = B 1 +B 2 



(30) 



where 



Si 



-V2 



7(£,*) 



Pa 



d 2 



1 



dz 2 4h 2 



[P-Pi 



P\P- 



(31) 



(32) 



in which 



1 (s, z ) = ^{p pJJ^-J^- 

3 L \GqGi EqE\ 



A 

T 



p2 
EqGq 



2 



1 

G~i 



P{ 
E\G\ 



2 

Go 



1 



(33) 

The diagonal components in the Hamiltonian (22) are 
composed of the terms proportional to Pg and P 2 , result- 
ing from the first iteration step (marked I in the previ- 
ous sections). They contain also the terms proportional 
to PqQ 2 and P 2 Q 2 , resulting from the third iteration 
step (marked III). The diagonal terms neither raise nor 
lower the harmonic oscillator functions. The nondiagonal 
terms B\ and B 2 are proportional to Q and they come 
from the second iteration step (marked II). They result 
from the bulk inversion asymmetry of the crystal. The 
operator B\ raises the harmonic oscillator function \n > 
to \n + 1 >, whereas B 2 lowers it from \n > to \n — 3 >. 
In consequence, solving the eigenenergy problem imposed 
by Eq. (22) requires again the Evtuhov procedure. How- 
ever, now we have to look for solutions in the form 



the growth direction z. When applying the functions (34) 
we will limit them to the minimal coupling scheme. It 
can be seen from Eq. (22) that the spin-up state In, + > 



> 

CD 



3 r 



G ^0.67 Al o. 33 AS 



GaAs 



G 3.0.67^1 0.33 AS 



, G, 



0.5- 



o - 



E, 



1.5 



-2 



Go 




E ^ 



50 







z(A) 



100 150 



FIG. 4: Band structure of a rectangular GaAs/Gao.67 AI0.33AS 
quantum well within 5LM along the growth direction z. The 
calculated energy and wavefunction of the ground subband is 
indicated. 



described by A + interacts via the offdiagonal elements 
Bi and B 2 with the spin-down states. In view of the 
above considerations it is clear that \n,+ > interacts via 
B 2 with \n + 3, — > and via B\ with \n — 1, — >. Af- 
ter performing the operations on the harmonic oscillator 
functions we obtain the following eigenenergy problem in 
the minimal coupling scheme 



Xn{z) 

Xn+-i{z) 1=0. 
Xn-l{z) 

(35) 

This set of coupled differential equations gives three 
energies, of which we are interested in £{n, +). A similar 
reasoning leads to the following set of differential equa- 
tions containing the \n, — > state 




(34) 



where 1=1, 2, the functions \m >= A m exp(ik x x)$ m (y) 
are the same as for the bulk, but Xm(z) are as yet un- 
known envelope functions describing the motion along 



A + 



B. 



t 

2,n 



£ B 2 . n 
A+ +1 -£ Bi, n 



B 



t 

1,71 



A,. 



Xn-3(z) 
Xn+l(z) I = 0. 
Xn(z) 

(36) 



Here we are interested in £ (n, — ). 

In the above matrices we use the following notation 
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j-i- • « ft 2 d 1 9 fte_B , . . iibB . 



+EQ ~t\2^Tz [EP " ((n + l)Kl + nK2) + EPl ((n + 1)K3 + nKi)] d~z + 



-hu>° c [E Po ((n + l) 2 if 5 + n(n - 1)K 6 + (n + l)(n + 2)Jf 7 ) + ^ ((« + l) 2 ^s + n(n - l)tf 9 + (n + l)(n + 2)X 10 )] } 

(37) 

A « = ^ - Tai^) ai + ^7) {n + 1/2) " z ) + 



M f h 9 

6 \ 2mo dz 



[E Po (nK, + (n + 1)K 2 ) + E Pl (nK 3 + (n + 1)K 4 )} 



dz 



~hw° c [E Po (n 2 K 5 + ( n + l)(n+ 2)K 6 + n(n - 1)K 7 ) + E Pl (n 2 K 8 + (n + l)(n + 2)X 9 + n(n - 1)K W )] } (38) 



B 1 ,„= 7 (£,z)^(^V^TT^ + 



(n + 1 
2L 2 



2,n 



V2L 3 



v/(n + 3)(n + 2)(n + l) 



(39) 



(40) 



in which K\, ...Kio are again given by Eq. (29). 

If one is interested in the lowest spin states |0, ± >, 
sets (35) and (36) do not contain the n — 1 and n — 3 
components and the eigenenergy problems reduce to the 
following two sets of coupled differential equations. For 
the 10, + > state 



/ Aq — £ _B 2 .n 

while for the |0, — > state 



Ai — £ Bi in 



Xo(z) 



= 



A -£ 



Xi(z)\_ 
Xo(z) 



(41) 



(42) 



In order to determine the cyclotron mass for the lowest 
levels we need to know in addition the energies of 1 1 , + > 
and |1,— > states. As follows from Eq. (36), for the 
|1, — > state the corresponding differential set still con- 
tains two equations, while for the |1, + > state we have 
in principle three coupled equations, see Eq. (35). How- 
ever, we found that the addition of the third state |4, — > 
(see Eq. (35)) does not really change the energy of the 



|1,+ > state, so that two equations suffice to determine 
its energy. 

Once the electron is in a QW of a finite height, its 
total g*-factor is given not only by its properties in the 
well but also by those in the barriers. The situation for 
a rectangular GaAs/Gao.67Alo.33As QW is shown in Fig. 
4. Since the electron wave function penetrates into the 
barriers in which g* is different from those in the well, 
we should average over the two regions. Seemingly, one 
should simply use the band edge values of g* in GaAs 
and Gao.67Alo.33As and take an appropriate average. In 
fact, this was done in Ref. 3. We argue, however, that 
this is incorrect. The reason is that, as follows from Fig. 
4, the electron energy is much lower than the barrier's 
band edge Vb- This effect is automatically included in 
our formalism if we use for both the well and the bar- 
riers the appropriate functions Ei and as given by 
Eqs. (27) and (28), respectively. When considering the 
two regions we have to change not only the appropriate 
values of energy gaps and spin-orbit energies but also the 
potential V(z). At the flat bottom of the well shown in 
Fig. 4, there is V(z)=0. On the other hand, in the bar- 
rier regions there is V(z) = Vb, where Vb is the offset 
value for the conduction band between the two materials. 
Clearly, the energy £ is the same in the two regions. Only 
if Vb were zero would we deal with "straight" values of 
g* in both materials. 

It is instructive to consider a numerical example. Tak- 
ing for Gao.67Alo.33As the parameters given in Table 1 we 
obtain the band edge value of g$ = 0.5378, in agreement 
with the experimental findings [13]. On the other hand, 
when we compute the value of g* in the barrier using the 
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same parameters but accounting for Vb = 0.264meV and 
taking to a good approximation £ w 0, we obtain g* — 
0.1638, which is a considerably lower value. 

In fact, in our formalism we do not employ the aver- 
aging procedure. When calculating the total electron g* 
value we use in the expressions for g}(£, z), m*j(£, z) and 
Ki, as given by Eqs (37), (38) and (29), the actual func- 
tions Ei and Gi for a given region. This includes V(z) 
and the local values of band parameters. Thus the ac- 
tual values of gf, m*j and Ki enter both the differential 
equations defined by Eq. (22), as well as the boundary 
conditions at the interfaces. The problem of the bound- 
ary conditions was intensely discussed in the early days 
of heterostructure investigations, see [4]. It will suffice to 
say here that our Eqs. (37) and (38) are written in the 
hermitian form, so that the boundary conditions are ob- 
tained in the form of continuity of the wave function and 
of electric current across the interfaces. The eigenenergy 
calculation is done for each spin separately so that, as fol- 
lows from Eq. (22), also the boundary conditions depend 
on the spin. Since the effective equations for the energies 
contain only p 2 , terms, the boundary conditions are not 
more complicated than those for a z-dependent effective 
mass. Once £q and £q are computed, the g*-value is 
determined using the definition (15). 

All the above remarks apply also to the computation 
of the cyclotron mass. The mass value in barriers also 
contains the offset Vb- In the above numerical example 
for Gao.67Alo.33As, instead of the band-edge value itiq = 
0.0875mo we deal with the value m* = 0.0748mo in the 
barrier. We emphasize that this feature is of importance 
even in the simplest calculation of the energy and of the 
effective mass in a QW made of two materials. The above 
problem is recognized in the work of Bastard et al [19, 
20] as well as in our earlier papers [21, 22]. 



IV. RESULTS AND DISCUSSION 

Figure 5 shows our calculation of the electron spin g*- 
factor in GaAs/Gao.67Alo.33 As rectangular QWs together 
with available experimental data. This is the main re- 
sult of our work. It can be seen that the complete five 
- level P p model gives an excellent description of the 
g*-value as a function of the well width between 3nm 
and 12nm. In particular, the theory reproduces the ex- 
perimental value of g* — for the well width of 6 nm, 
as observed by Le Jeune et al [23]. In the discrepancy 
of the experimental points around d w 4.5 nm, priority 
should be given to Ref. [23] (the full circle) since the 
Raman data of Ref. [24] (the empty circle) do not deter- 
mine the conduction g*-value directly. Our calculation is 
carried out for the growth direction z||[001] with an ex- 
ternal magnetic field parallel to this direction. In order 
to reach the good description of experiments we had to 
use the full five-level model of the band structure, as dis- 
cussed above. The theory agrees perfectly well with the 
experimental <?*-value even for d sa 30 A. One could find 
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FIG. 5: Spin g*-factor of conduction electrons in rectangular 
GaAs/Gao.67 Alo.33 As quantum wells versus well width d. The 
solid line - theory based on 5LM, the dashed line - theory 
based on 3LM, full circles - experimental data of Le Jeune et 
al [23], empty circles - experimental data of Sapega et al [24]. 



it surprising that the P-p theory works so well on the 
scale of few interatomic distances in GaAs, but this result 
is in agreement with a more general experience that the 
P-p theory works better than it should. It can be seen 
that the total g*-value for narrow wells is smaller than 
the value for pure GaAs, as shown in Fig. 5. This is 
due to the above mentioned effect of the offset Vb , which 
lowers the effective g* in the barriers, as compared to the 
band-edge <?*-value in Gao.67Alo.33As. 

The dashed line in Fig. 5 shows the calculation of g*- 
factor, in which we used the three-level P-p model. This 
model does not involve the Q matrix element, so that 
it includes neither BIA splitting nor the Q 2 terms. It 
can be seen that 3LM does not properly describe the ex- 
perimental data in GaAs/Ga . 67 Alo.33 As quantum wells. 
This result agrees with the general conclusion that the 
three-level model does not describe correctly the conduc- 
tion band in GaAs - type materials, see [1, 2]. The values 
of g* shown in Fig. 5 were computed for the magnetic 
field intensity B = 2 T. It is natural that in a theory 
for the nonparabolic energy band, when the energies are 
nonlinear functions of a magnetic field, the g*-value may 
depend somewhat on B. We checked this dependence 
and found that diminishing the field until B w 0.25 T 
the g*-value varies only weakly with B. However, as B 
tends to zero the g*-value ceases to be a useful quantity 
because of the spin splitting due to BIA. It is then more 
practical to use the energy splitting for the two spins. 

As to the previous description of g*-value in QW based 
on 3LM [3], it is incorrect on two accounts. First, as al- 
ready mentioned in the discussion of Fig. 2, it predicts a 
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FIG. 6: Cyclotron mass of conduction electrons in rectangular 
GaAs/Gao.74Alo.26As quantum wells versus well width d. The 
solid line - theory based on 5LM, the dashed line - theory, 
based on 3LM, experimental data are after Huant et al [25] 
for samples with slightly different chemical composition. Full 
circles x—0.26, full square x=0.25, full triangle x=0.23. 



much too strong increase of g* in GaAs with decreasing 
well width d. Second, when averaging the spin g*-factors 
in GaAs well and Gao.67Alo.33As barriers, it takes the 
band-edge values in both materials, neglecting the influ- 
ence of the band offset on the <?*-value in barriers. Both 
above shortcomings lead to an increase of the total g*. 
We conclude that the agreement of the theory [3] with 
the experiment on GaAsZGao.67Alo.33As QW, as quoted 
by Le Jeune et al [23], is fortuitous. 

In Fig. 6 we plot the cyclotron mass of electrons in 
GaAs/Gao.74Al .26As rectangular QWs as a function of 
the well width according to the three-level and five-level 
models. It can be seen that 3LM gives distinctly lower 
masses, although we adjust the far-band contribution C 
to get the same band-edge value of = 0.0665 mo. The 
five-level model gives slightly different masses for spin-up 
and spin-down cyclotron transitions. 

Figure 6 shows the comparison of our theory with the 
cyclotron resonance data of Huant et al [25] obtained on 
GaAs/Gai_ x Al x As QWs. For the GaAs/Gao.7 4 Al .26As 
alloy we use the band parameters indicated in Table 1. 
For the band edge in GaAs we take the perfectly accept- 
able value of ttiq — 0.0665 mo, accounting for the fact 
that the data of Ref. [25] were taken at the tempera- 
tures 20-30 K. Our theory describes the data very well 
with the exception of the point at d w 20 nm, but this 
point does not follow well the overall mass behavior. The 
calculation shown in Fig. 6 was performed for B = 10 T, 
which corresponds to the experimental conditions of Ref. 
[25]. Comparing the data with the theory we conclude 



that the three-level model can not correctly describe the 
experimental data. 

The previous theoretical treatment of the effective 
masses in GaAs/Gai_ x Al x As QWs by Ekenberg [4], 
which was carried out for B = 0, makes a distinction 
between the masses parallel and perpendicular to the 
growth direction. We do not introduce the parallel mass 
since it does not seem to have a clear physical meaning. 
Also, in nonparabolic energy bands one can not separate 
the electron motion in different directions since an in- 
crease of the mass due to motion in one direction will 
affect the motion in other directions. The theory of Ref. 
[4] is based on an expansion of the energy in powers of 
momentum [26]. The coefficients in such an expansion 
contain the band-edge energies in their denominators and 
the interband matrix elements of momentum in their nu- 
merators. In Ref. [4] the potential V(z) of the well is 
added later, so that this procedure can be qualified as 
semiclassical. On the other hand, in our treatment we 
introduce the potential V(z) and the energy £ from the 
beginning and use the iteration procedure in solving the 
resulting 14 equations by substitution. This procedure 
is exact for the matrix element Q — and includes Q 
to second order. The resulting effective equations have 
the coefficients which not only contain the band-edge en- 
ergies, but also involve in their denominators V{z) and 
£, cf. Eqs. (27), (28), (29). We showed above that 
these quantities are not negligible in GaAs-type materi- 
als and they would be even more important in narrow- 
gap materials. The work [4] does not include the bulk 
inversion asymmetry effects, which are not very impor- 
tant for the cyclotron mass but become important for 
the spin g*-value, see our Fig. 3. Our theoretical masses 
for GaAs/Gao.7Alo.3As QW (not shown) are somewhat 
higher then those presented in Fig. 4 of Ref. [4] . It is not 
clear whether the differences result from the approxima- 
tions made in Ref. [4] or from somewhat different band 
parameters used in the two calculations. 

The approach presented above can be used equally well 
for quantum wells described by an arbitrary potential 
V(z). This amounts to solving differential equations (22), 
(37), (38), in which V(z) explicitly appears. However, if 
a well (and a corresponding potential) does not possess 
the inversion symmetry, there appears a spin splitting 
due to the Bychkov-Rashba mechanism [6] . Technically, 
this appears as a result of noncommutativity of p z with 
V(z) when one solves the P-p equations by substitution. 
The Bychkov-Rashba splitting has been treated in many 
papers, see the review [27]. 

The present work confirms our previous description of 
the spin c/*-factor in parabolic GaAs/Gai-xAl^As quan- 
tum wells [28, 29]. We conclude that the five-level P-p 
model adequately describes this important system. 
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V. SUMMARY 

We describe the spin and cyclotron energies of elec- 
trons in GaAs/Ga^j-ALrAs quantum wells in the pres- 
ence of an external magnetic field parallel to the growth 
direction [001]. Our approach is based on the five-level 
Pp model of band structure for GaAs-type materials. 
Inadequacy of the previous theory used in the literature 
for the spin g*-factor in heterostructures is indicated. 
We solve 14 coupled differential equations resulting from 
the P-p formulation by three iteration steps. The suffi- 
ciency of this procedure is tested on bulk GaAs energies 
for different values of k\ corresponding to real quantum 
wells. Influence of the bulk inversion asymmetry present 



in III-V compounds on the electron g*-value is empha- 
sized. Our theory gives an excellent description of the 
existing experimental data on the spin g*-factor and the 
cyclotron mass of electrons in GaAs/Gai_ x Al x As rect- 
angular quantum wells for different well widths. 
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